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The temperature dependence of the resistance -R(T) of 
ultrathin quench-condensed films of Ag, Bi, Pb and Pd has 
been investigated. In the most resistive films films, R(T) — 
Ro exp (To/T) x , where x = 0.75 ± 0.05. Surprisingly, the ex- 
ponent x was found to be constant for a wide range of Ro 
and To in all four materials, possibly implying a consistent 
underlying conduction mechanism. The results are discussed 
in terms of several different models of hopping conduction. 



I. INTRODUCTION 

In highly disordered materials, electrical conduction 
occurs by the hopping of electrons between localized sites. 
This results in a thermally activated electrical resistance 
of the form: 



R(T) = 1. . .(^ 



(1) 



where T is temperature, and Rq, Tq and x are constants 
which depend on the disorder, the details of the inter- 
actions and the dimensionality of the system. Simple 
activated hopping over a constant barrier results in the 
Arrhenius form with x = 1. For noninteracting electrons, 
when the average hopping distance depends on temper- 
ature due to the compromise between hopping to sites 
which are close in energy, but farther away, Mott variable 
range hopping 0] is expected, with x = l/(d+l), where d 
is the dimension. Efros and Shklovskii (ES) showed that 
including Coulomb interactions between electrons results 
in a soft gap in the density of states at the Fermi energy, 
which changes the variable range hopping exponent to 
x — 1/2 in all dimensions 0. 

Hopping conduction has been investigated in a wide 
variety of materials, such as doped semiconductors HQ, 
semiconducting heterostructures j^j, amorphous metals 
fj|-|l0f , magnetic materials jll| and superconductors jl2) . 
Both the Mott and the ES forms of variable range hop- 
ping have been observed, as well as a crossover between 
the two regimes ||||. It should be emphasized, how- 
ever, that it is often hard to distinguish between Mott 
and ES hopping, particularly in experiments in which the 
resistance changes only by one or two orders of magni- 
tude. The unambiguous identification of the Mott or ES 
hopping can be further complicated by factors which are 
usually neglected, such as the granularity of the system, 



possible temperature dependence of Ro, or correlations 
between electron hops. 

While investigating the transport properties of ultra- 
thin quench-condensed films over the course of many 
years, we have often found that the resistance of the 
thinnest films was thermally activated with x ~ 0.75. 
A similar hopping exponent has been reported by other 
authors ||,|H],|l3],[l4| , but has rarely been discussed per 
se. Since there is no theory which predicts this value of 
the hopping exponent, its origin has been left an open 
question. Here we report a detailed study of the tem- 
perature dependence of the resistance in very disordered 
films of four different materials: Ag, Bi, Pb and Pd. The 
films were grown in separate runs over several years, in 
different cryostats and on different substrates, yet they 
all show the same thermally activated resistance with an 
almost identical exponent. A careful analysis of the data 
points to a new conduction mechanism in this regime, or 
perhaps calls for a modification of the conventional pic- 
ture. We compare our results with those of other experi- 
ments and available theoretical calculations, and suggest 
that the model that may possibly explain the anomalous 
hopping exponent is the collective variable range hop- 
ping model of Fisher et al. jLj|. Developed to describe 
vortices in superconductors, this approach has not been 
considered before in the context of charge transport in 
disordered electronic systems. 

In Section II we survey the various models that have 
been considered in the discussion of transport in disor- 
dered films. We also exhibit the mapping of the model 
for collective vortex hopping onto the problem of charge 
transport in disordered systems and estimate the value 
of the hopping exponent. Experimental details of film 
growth and resistance measurements are given in Sec. 
III. In Sec. IV, we analyze the temperature dependence 
of the film resistance using several different methods to 
show that the exponent obtained is really a property of 
the system, rather than a consequence of an improper 
fit. The results are discussed and compared with other 
experimental and theoretical work in Sec. V. 



II. SURVEY OF HOPPING MODELS 

In recent years, extensions of the basic variable range 
hopping model to include percolation effects and cor- 
relations between electron hops have been developed. 
Deutscher et al. fl6[| proposed a hopping mechanism 
which leads to a thermally activated resistance with x 
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close to 1/2 without considering Coulomb interactions. 
The mechanism was based on the superlocalization prop- 
erty of wavefunctions on incipient percolation clusters 
|p7| , and may be relevant for atomically disordered sys- 
tems as well as for granular percolative structures. The 
detailed microstructurc of ultrathin quench-condensed 
films is still not known, and although these films are 
usually considered to be homogeneous, they may actu- 
ally contain small grains or clusters. It is then possible 
that the electrons are restricted to move on a sublattice 
which is fractal over some range of length scales, and 
that their wavefunctions decay faster then exponentially 
with distance. Based on this assumption, a hopping con- 
ductivity law has been derived [|l6| near the percolation 
threshold, which has the form of Eq |l| with x = 3/7. 
This is experimentally almost indistinguishable from the 
Efros-Shklovskii law with x — 1/2 if only the temper- 
ature dependence of the resistance is studied, but can 
be identified through the behavior of the parameter Tq 
and the nature of the crossover to the conventional Mott 
regime fl6|| . 

Generally speaking, there is no reason to assume that 
the prefactor Rq in Eq. Jh is independent of tempera- 
ture. Van Keuls et al. g studied the resistivity in a 
gated GaAs / Al x Gai- x As heterostructures as a function 
of temperature, electron density and magnetic field. As- 
suming the prefactor of the form: 



Ro = bT n 



(2) 



where b and m are constants, these workers fit Eq. |l] to 
their data with x = 1/3 in low magnetic fields and with 
x = 1/2 in high magnetic fields. The same crossover 
was observed as a function of electron density and tem- 
perature, and it scaled with the separation between the 
electron layer and the nearby screening gate, as predicted 
by Aleiner and Shklovskii |Tj| . In addition to introducing 
a temperature dependent prefactor, this experiment also 
raised the issue of the importance of correlations between 
the electron hops. 

In an electron glass, where the screening length is long 
p9[ and the interactions long-ranged, electron hopping 
may be correlated [^0| . Excitations can leave the system 
far from equilibrium pl| and relaxation occurs through 
the rearrangement of charge. The energy of a single elec- 
tron hop may then be significantly reduced by the motion 
of the surrounding charges. At sufficiently low temper- 
atures, such collective hopping might be the dominant 
conduction mechanism. 

In the analysis of their data, Van Keuls et al. || as- 
sumed that the number of configurations of occupied 
states reached by the correlated hopping of a number 
of electrons is proportional to the single-particle density 
of states. In that case, the qualitative behavior of the 
resistance remains unchanged, and the effects of corre- 
lations enter through the constants which determine the 
parameter To in different regimes. 

Yet another issue which can be relevant in extremely 



thin films is the possibility that the electrons might in- 
teract logarithmically rather than as 1/r. As shown by 
Keldysh [^2| , the range of the logarithmic interaction is 
given by: 



1 



(3) 



where n and the dielectric constants of the film and 

the substrate, respectively, and d is the film thickness. 
One possible consequence of the logarithmic electron- 
electron interactions is that the Coulomb gap in the den- 
sity of states (linear in 2D) may change to an exponential 
form p3fl . This leads to a modified variable range hop- 
ping law with a temperature dependent exponent. 

Alternatively, the behavior of the logarithmically in- 
teracting electrons might be similar to that of vortices in 
2D superconductors, which are known to interact loga- 
rithmically. Collective variable range hopping of vortices 
in disordered thin-film superconductors was studied by 
Fisher et al. Jis] , Including the effect of correlations, 
these authors found that multivortex hopping results in 
a lower energy than single vortex hopping. They also 
suggested that such multiparticle hopping might domi- 
nate single particle hopping even in the case of inverse 
power law interactions. 

This approach may then be mapped onto a disordered 
two-dimensional system of charges. Following the argu- 
ments of Fisher et al. [ ]l5| |, the energy U{r) of a multi- 
particle excitation of length r, can be estimated to be: 



U(r) « K{-) 1 / 2 



(4) 



where I is the distance between charges, and K is the 
bare single-particle excitation energy. The latter is equal 
to e 2 / na, where a is the localization length and k is the 
dielectric constant. The simultaneous hopping of many 
charges may then result in a lower energy than the hop 
of a single charge. 

The electrical resistance is a product of the probability 
for an electron to tunnel a distance r, T r , and the prob- 
ability for an excitation with the energy U(r) to occur, 



Rot l/T r T L 



(5) 



A lower bound for the multihop rate can be estimated 
as follows: if all of the (r/l) 2 electrons in the r by r region 
hop a distance comparable to the spacing I, then the rate 
should be proportional to the single-hop rate , exp(—l/a), 
raised to the power of the number of electrons, resulting 



-(l/a)(r/lf 



(0) 



The probability for an excitation of energy U (r) is pro- 
portional to exp(—U(r)/T), or using Eq. W: 



oc e 



-(K/TKl/r) 1 ' 2 



(7) 
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The minimum resistance is obtained when the hopping 
distance r^ op is: 

rl P « aUf ) 4/5 (8) 

where 

T = K (//a) 1 / 4 (9) 

Substituting Eqs. || through || back into Eq. [j] results 
in a hopping form such as that of Eq. ^ with x = 4/5. 

In the other limit, the minimum number of electrons 
participating in a collective hop could be taken as (r/d), 
which leads to x = 2/3. The collective variable range 
hopping mechanism may therefore result in a resistance 
of the form of Eq. [l], where the range of exponents is 
2/3 < x < 4/5, depending on the fraction of electrons 
participating in the process. 



a linear function of the applied current. When measur- 
ing very resistive films (10 4 — 10 s f2), because of the long 
time constants of the circuit, it was necessary to mon- 
itor the voltage as a function of time after the current 
was changed, and allow adequate time for the voltage to 
stabilize. To avoid the voltage offset errors due to ther- 
mal EMFs, both polarities of the current were used to 
determine the resistance. 

The resistance of a series of Ag films was also studied 
in a magnetic field. Magnetic fields of up to 20kG (12 
kG) were applied in direction parallel (perpendicular) to 
the plane of the substrate using a split-coil superconduct- 
ing magnet. In a regime where the anomalous hopping 
exponent is observed, the resistance was found to be in- 
dependent of magnetic field. 

IV. RESULTS AND ANALYSIS 



III. EXPERIMENTAL METHODS 

The temperature dependence of the resistance has been 
studied in ultrathin quench-condensed films of Ag, Bi, Pd 
H and Pb (24|j2f|. The films were deposited on liquid 
helium cooled substrates and resistance measurements 
were performed in situ at temperatures down to 0.15/V 
Ultra-high vacuum conditions and temperatures below 
20K were sustained throughout each run, in order to 
avoid contamination or crystallization. The substrates 
were glazed alumina (for Bi and Pd films) or SrTiO^ 
(100) (for Ag and Pb films). The SrTi0 3 (100) sub- 
strates were 0.75mm thick and had a lOOnm thick Au 
gate on the back. Such a gate can be used to study the 
response of the film to a perpendicular electric field and 
was used to establ ish glass-like behavior in the most re- 
sistive films 
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Films were deposited in thickness increments between 
0.05 and 0.5 A on top of a thin germanium layer (5— 10 A). 
(The Pd films were the exception to this as they were 
deposited directly onto glazed alumina substrates where 
they became connected at monolayer coverage.) Films 
grown on amorphous Ge are believed to be homogeneous, 
since they become connected at an average thickness of 
about one monolayer [p6[ . The thicknesses of the films 
studied ranged from bA up to 15A. These nominal values 
of the film thickness are calculated from the deposition 
rate, which was measured using a quartz crystal mon- 
itor placed in the vicinity of the substrate. The first 
low-temperature scanning tunneling microscopy studies 
of the morphology of films grown in a similar manner in- 
dicate that the thinnest films may indeed be homogenous, 
while the thicker ones contain small clusters (27j . 

Resistance measurements were carried out using a 
standard dc four-probe technique. Very low bias currents 
(< lOnA) were used to avoid Joule heating of the sample 
and to make sure that the voltage across the sample was 



The temperature dependence of the sheet resistance 
(resistance per square) for series of Ag, Bi, Pb and Pd 
films is shown in Fig. EL The logarithm of the sheet resis- 
tance, plotted as a function of T -0 75 follows a straight 
line for each film, indicating that the resistance is ther- 
mally activated with x « 0.75. Using other values for x, 
such as 1, 1/2 or 1/3 yielded considerably larger devia- 
tions from the data. 

Since the prefactor in Eq. [j] is generally expected to be 
temperature-dependent (m. ^ in Eq. ^|), we attempted 
to fit the data using different combinations of m and x. 
As shown in Fig. ^, using values of m greater than zero 
actually increased the error of the fit. The maximum 
deviation in the fit of the combinations of Eqs. 1 and 2 
to the data became much larger than the noise in R as 
m was increased. Furthermore if values of m were chosen 
to force either Mott or ES hopping exponents of x = 1/3 
or 1/2, respectively, the quality of fits as measured by 
chi squared would be significantly worse than that with 
m = 0, in contrast with the findings of Van Keuls et al. 

Assuming that the hopping exponent is i k 0.75, the 
activation energy To can be extracted from the fit to Eq. 
|l|. The values of the parameter To for different films of 
all four materials are shown in Fig. ^| as a function of 
Ri4K, which is the sheet resistance measured at 1AK. 
This quantity is inversely proportional to the film thick- 
ness, so by using Ruk instead of the thickness, one can 
avoid systematic errors in the nominal thicknesses of the 
films of different materials. It is apparent in Fig. || that 
To changes greatly as Ruk (and therefore also the film 
thickness) changes, ranging from around 100K for the 
thinnest films, to around 1QK for the thickest films. The 
same qualitative and quantitative behavior was found for 
all four materials. 

A more direct method of determining the exponent 
x (which is exact under the condition m << (Tq/T) ) 
has been developed by Zabrodskii and Zinov'eva [EJ. 
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The method is based on defining the function w = 
— d(log R) / d(logT) . If R is given by Eq. [I], then logw cx 
— xlogT. By plotting logw as a function of logT, x can 
be easily extracted from the slope of the resulting straight 
line. The benefit of the Zabrodskii-Zinov'eva approach 
is the simplicity of fitting a line rather than a compli- 
cated function with up to four adjustable parameters(6, 
m, Tq and x). Once it has been determined that m = 0, 
this method eliminates the danger of finding a local min- 
imum instead of the best fit. The results of determining 
x this way are shown in Fig. i| . For very resistive films, 
plotting log w vs. log T indeed yielded straight lines. The 
values of x varied slightly between different materials, 
from x w 0.7 for Ag to x « 0.8 for Bi. Remarkably, the 
value of x did not change between different films in the 
same series over a significant range of sheet resistances, 
as shown in Fig. [5[ 

Even though Ruk and To change from film to film, as 
more material is added to increase the average thickness 
and decrease the sheet resistance of the film, x stays con- 
stant over three orders of magnitude in R\^k- For thicker 
films (smaller R\^k)i % drops rather abruptly to a value 
between 1/3 and 1/2. In this regime, the Zabrodskii- 
Zinov'eva plots no longer produce straight lines, indicat- 
ing that the hopping exponent changes as a function of 
temperature. Further increase of the film thickness leads 
to another crossover to a weakly localized regime where 
the temperature dependence of the resistance is logarith- 
mic (not shown in Fig. H). 



V. DISCUSSION 

An activated temperature dependence of the resistance 
with an anomalous hopping exponent x ks 0.75 has been 
observed in disordered films of four different materials, 
grown on different substrates and measured in different 
cryostats. This strongly suggests that the exponent x f» 
0.75 is a general property of ultrathin films of metals in 
the very strongly localized regime. 

The same exponent has been reported by Adkins and 
Astrakharchik jl0| in ultrathin quench-condensed films 
of Bi with a Ge underlayer. In that experiment, the tem- 
perature dependence of the resistance changed to simply 
activated (with x = 1) when the Coulomb interaction 
was screened in the presence of a nearby metallic gate. 
This behavior was ascribed to the fixed range hopping 
of dipoles in screened films, but no details were given on 
the origin of a; « 0.75 in unscreened films. The authors 
suggest that the films may be in the crossover regime 
between the variable range hopping and the fixed range 
hopping regime. Such a crossover can occur when the 
optimal hopping distance r^op becomes comparable with 
the localization length £. 

In our experiment, x stays constant over several or- 
ders of magnitude in sheet resistance, and then drops 
abruptly as the sheet resistance decreases further. If our 



films were merely at the crossover between x ~ 1/2 (or 
1/3) and x « 1, the change in x would be expected to 
be gradual. The observed constancy of x implies that 
a consistent mechanism may be governing the conduc- 
tion in this regime, a mechanism different from Mott or 
Efros-Shklovskii variable range hopping which is usually 
observed in less resistive films. 

Furthermore, it was not possible to obtain a satisfac- 
tory fit to the data using a temperature-dependent pref- 
actor, as in the work of Van Keuls et al. It is interesting 
that these workers obtain x « 0.75 in all magnetic fields 
if m is taken to be zero. However, the activation energies 
obtained from such fits are reported to be unacceptably 
small || . 

There are several other mechanisms which may be rele- 
vant in a very disordered 2D system. For example, Dai et 
al. H observed an exponent x ~1 in Si:B, which changed 
to x ~l/2 when a magnetic field was applied. They sug- 
gested that the x ~1 was due to the exchange interaction 
between the electron spins, which is destroyed in a mag- 
netic field. It must be noted that the Ag did not show 
any magnetoresistance up to the highest field available, 
20 kG, so the exchange interaction is most likely not the 
origin of the anomalous hopping exponent observed in 
these films. 

If we allow a possibility that our films are granular 
on a very small scale (which we cannot unequivocally 
rule out), than we must consider the superlocalization 
mechanism of Deutscher et al. (ll| as a possible candi- 
date to explain our data. Without Coulomb interactions, 
this model predicts x « 0.43, which obviously cannot ac- 
count for our results. Including the Coulomb interactions 
may lead to a higher exponent, as proposed by van der 
Putten et al fllij . These authors studied the hopping 
conductivity of percolating carbon-black-polymer com- 
pounds and found x m 0.66, which is much closer to our 
result, although still too low. They interpret their re- 
sults as evidence of superlocalization on a fractal network 
with Coulomb dominated hopping. The activation ener- 
gies were found to be independent of the electron con- 
centration, as predicted by Deutscher et al. In contrast, 
the activation energies found in our experiment depend 
strongly on the film thickness, which is closely related to 
the electron concentration. Furthermore, if the Coulomb 
interactions were screened, one might expect the expo- 
nent to decrease towards 0.43, rather than to increase 
towards 1, as observed in screened Bi films by Adkins 
and Astrakharchik Jl(J . 

Another possibility is that the anomalous exponent is 
a consequence of the exponential gap in the density of 
states, which can arise if the electrons interact logarith- 
mically p3| . In that case, the hopping exponent would 
be something close to, but smaller than 1 at higher tem- 
peratures, and then cross over smoothly to 1/2 at low 
temperatures. Forcing a fit of Eq. [I] to the data would 
result in an exponent which changes continuously with 
temperature. In the less resistive regime where we ob- 
serve a temperature-dependent exponent, a closer inspec- 
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tion shows the opposite trend: the exponent is close to 
1/2 at higher temperatures, and increases with decreas- 
ing temperature. On the other hand, we cannot rule out 
the possibility that we might observe a smooth crossover 
to 1/2 in the most resistive films if we could measure at 
much lower temperatures. 

Finally, we consider the collective variable range hop- 
ping mechanism, proposed in the context of vortices in 
disordered superconductors by Fisher et al. |lq ]. The 
mapping of this model onto a 2D electron system may 
actually be exact, if the electron-electron interactions are 
logarithmic over relevant length scales, but the authors 
suggest that collective hopping may dominate over the 
single-particle hopping even in the case of a conventional 
Coulomb interaction. The range of exponents predicted 
by the collective hopping model is 2/3 < x < 4/5, de- 
pending on the ratio of electrons which participate in 
the process. The exponent found in our experiment, as 
well as the exponents found by other groups |l^,||,|l^,|lj| , 
are well within that range. The activation energies are 
expected to depend on the localization length, which in 
turn depends on the film thickness, as shown in Fig. ^. 

In conclusion, we have addressed the issue of the 
anomalous hopping exponent x w 0.75 observed in ultra- 
thin films of metals and related 2D systems. We argue 
that this hopping exponent is a general property of very 
strongly disordered systems, rather than a result of an 
improper fit or a signature of some sort of a crossover 
behavior. The usual models of hopping conduction do 
not explain this result. Our data can be explained by 
a collective variable range hopping mechanism, but our 
work by no means provides a proof of such a mecha- 
nism. Future experimental and theoretical studies will 
be needed to shed more light on this matter. 
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FIG. 1. Sheet resistance as a function of T~ ' 75 for a series 
of Ag, Bi, Pb and Pd films. Each curve represents a different 
thickness, and the film thickness increases from top to bottom: 
[1] N. F. Mott, J. Non-Crystal. Solids 1, 1 (1969). a) Ag 5 _ 6 . 6 A, b) Bi 9-10A, c) Pb 4.9-5.2A, and d) Pd 10-13A. 

[2] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 
[3] P. Dai, Y. Zhang, and M. P. Sarachik, Phys. Rev. Lett. 
69, 1804 (1992). 
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FIG. 2. Best fits for x (full circles) for increasing exponent 
m in the resistance prefactor T m in the combination of Eqs. 
1 and 2 for a 7. OA thick Ag film. Chi squared (open circles) 
is a minimum at m = 0. The deviation of the fit from the 
data quickly becomes much larger than the noise in R as m 
is increased. 

FIG. 3. The activation energy To as a function Riak for 
different films of all four materials: Ag (crosses), Bi (squares), 
Pb (diamonds) and Pd (circles). 

FIG. 4. Zabrodskii plots for a series of Ag, Bi, Pb and Pd 
films. Each curve represents a different thickness, and the 
film thickness increases from top to bottom: a) Ag 5-6. 6A, 
b) Bi 9-10A, c) Pb 4.9-5.2A, and d) Pd 10-13A. The slopes 
determine the exponent x. 

FIG. 5. The hopping exponent i as a function of the 
sheet resistance Ri4K for all four materials: Ag (crosses), Bi 
(squares), Pb (diamonds) and Pd (circles). 
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